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We present an exact solution to an interface model representing the dynamics of a 
domain wall in a two-phase Ising system. The model is microscopically motivated, 
yet we find that in the scaling regime our results are consistent with those obtained 
previously from a phenomenological, coarse-grained Langevin approach. 
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The study of the magnetisation profile between coexistent phases in mod- 
els of subcritical uniaxial ferromagnets (and their analogues) has generated 
important and unexpected theoretical advances. A typical example is the 
existence of spatial fluctuations of the interface location so large that they di- 
verge with system size. Phenomenologically speaking, these fluctuations are 
generated by capillary waves described in a continuum theory of an interface 
which is a sharp dividing surface between oppositely magnetised phases, con- 
trolled by a somewhat arbitrary short range spatial cutoff [[]]]. Against this, 
one has to set the density functional theories going back to van der Waals 
and Maxwell [[| which do not give divergent spatial interface fluctuations. It 
is fortunate that both exact and rigorous results are available for the planar 
Ising model which resolve the conflict in favour of the large spatial fluctua- 
tions ||, but have yet to resolve the problem of local interface structure. The 
conclusions of the exact calculation were recaptured in a Helmholtz fluctua- 
tion theory for the phase separating surface [|J which brings in the concept 
of interfacial stiffness and which is valid in a spatially coarse-grained sense, 
much as in Widom scaling theory ||. These fluctuations are an essential 
ingredient in the statistical mechanics of a number of surface phase transi- 
tion phenomena 0. In this letter we focus on dynamical aspects of these 
fluctuations by studying an exactly solvable model on a finite lattice. 

First, we recall an equilibrium result which helps motivate the model 
itself and an earlier phenomenological treatment. Suppose in a zero-field 
planar Ising ferromagnetic model, the interface between coexistent phases is 
established and localised in laboratory-fixed axes by specifying the boundary 
spins as shown in Fig. 1. The straight line of length 2L connecting the spin- 
flip points is the Wulff shape for the interface. It is convenient to define 
coordinates parallel and perpendicular to this line, with origin at its centre. 
Let the thermodynanic limit of infinite strip length be taken first. Then 
denoting the magnetisation at (x,y) by m(x,y/L) we have for — 1 < a < 1 
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where m* is the spontaneous magnetisation, <fr(x) is the error function erf (x) 
and &($, T) = er($) + is the surface stiffness. (Here is the 

angle-dependent surface tension.) On the other hand, one has m{a,(3) = 
(m*sign( / 9)) for 5 < 1/2 (> 1/2). 

In this letter we study the dynamical case with a new exact result. In 
order to generalise to the dynamical case we need to define a model with 
dynamics appropriate to the physical system described above. We shall see 
that our model reproduces the equilibrium result (|l|) in the infinite time 
scaling limit. To motivate our model, note that as T — > 0, &($) > (strictly) 
provided $ ^ 0,±tt/2; for & = 0,±?r/2, 6(0, T) -> as T -> 0. This is a 
primitive example of facetting. At T = 0, the Peierls contours reduce to paths 
connecting (-L, 0) to (L, 2Ltani?) which either step to the right or upwards 
(tan-# > 0); all such paths are degenerate. They separate regions of opposite 
magnetisation; since T = 0, this magnetisation is of unit magnitude. In this 
case, the theorem is related to normal fluctuations in coin tossing. Let us 
now concentrate on the case tan$ = 1, with an initial sawtooth configuration 
as shown in Fig. 2. For later convenience we have rotated the system by 45 
degrees. 

At T = 0, any minimum energy path can be represented as a sequence 
S = {n-L+i, n-L+2, ■ ■ ■ of 2L binary numbers where nk — 1 if the 
k th segment of the interface steps upwards (in an angle of 45 degrees) and 
n.fe = if the steps goes downwards (see Fig. 2). One can think of as an 
occupation number which is related to the interface height hk in the rotated 
system by 2rik — l = hk — hk-i- The configuration at any time t is then given 
as a time-slice S(t) = {rik(t) : k = —L + 1, . . . ,L} and the dynamics are 
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given by specifying rules connecting S(t + 1) to S(t). 

We take particle-hole exchange on neighbouring lattice sites with prob- 
ability p, but leave pairs of particles or holes unchanged. This corresponds 
to single spin-flips at local extrema of the interface. The nearest neighbour 
interaction makes it advisable to update on alternating sublattices at sub- 
sequent time-steps. This updating scheme turns our model into a cellular 
automaton which is related to the six vertex model on a lattice oriented at 
45 degrees 0, |9|: take Uk = 1 (resp. =0) to be an upwards (resp. downwards) 
pointing arrow as in Fig. 3 where we show the ice-vertices, the correspond- 
ing dynamical events and their weights in the particle-hole picture. With 
< p < 1, the two row-to-row transfer matrices in the (1,1) direction gener- 
ate a stochastic time evolution satisfying local detailed balance. 

The transfer matrix can be written as 

L-l L 

T = T even T odd = J| ^ JJ ^ 
j=-L+l j=-L+l 

with Tj = 1 -pej, ej = ( ar j< J j + i + c r j&j + i + a 'j cr j + i ~ I)/ 2 , the a x > y ' z being the 
Pauli matrices. We choose — (1 — o"|)/2, a vacuum |0) with rz^ 1 0) = 0. 
The initial state is the L-particle state \v) = (|1, 0, 1, 0, . . .) + |0, 1, 0, 1, . . .))/2. 
The displacement of the interface at (x, t) is thus 

X 

h(x,t)= (l-2n fc (t)) (3) 

k=-L+l 

with zero mean by spin-reversal symmetry, but with 

w 2 ( Xl t) = ((h(x,t)) 2 } =4 ]T (LKn ; T>)-(x + L) 2 . (4) 

k,l=-L 

Here (L\ is the out state which sums over all eigenstates of the with 
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particle number N = L, giving each such state equal weight. Evidently 



(L\ = (L\)- l (0\(S + ) L (5) 
where S ± = J2 s t an d s k = (of ± ^)/2. Using the commutator 

^+) L , n fc ] =L a +(5 + ) L " 1 (6) 
and the SU(2) symmetry of T, (|5]) and (|]) together with (|4]) give 
u; 2 (a;,t)=4 £ E M 1 (k,l-p,q\t)M 2 (p,q)+2(x + L)-(x + L) 2 (7) 

k,l=—L+l p,q=-L+l 

where 

Mx(A;,/;p,g|t) = (0|4 S +T' S ;^|0) (8) 
is a matrix element in the 2-particle sector and 

M 2 (p,q) = (0\s; S +(S + ) L - 2 \v)/(L-2)\ (9) 

contains the information coming from the initial condition. Using (^) and 
the expression for \v) above, it follows that M 2 {p,q) — 1/2 for (p — q) even 
and =0 otherwise. Thus 

x L 

w 2 (x,t) = 2 J2 E M 1 (k,l;p,q\t) + 2(x + L)-(x + L) 2 (10) 

k,l=— L+l p,q=-L+X 

p-geven 

which is an enormous simplification of @). Expressions for higher moments 
can be obtained in a similar way. 

We study first the infinite-time limit, thus deriving the static magnetisa- 
tion m(x, y) = 1 — 2(Q(h(x) —y)) where 0(s) is the Heavyside step function. 
The (normalised) L-particle steady state of the system is \L)/B(2L,L) with 



the binomial coefficient B(m,n) = m\/(n\(m — n)\). A short calculation 
yields 



. . sign(y) ^ _. T L-x , . _,_ L + x , . . . 
™fay) = R/tr n E B ( L -x,^ — -k)B(L + x,— — + k) . (11) 

D\IL,L) k= _ y+l I I 

In the scaling limit a; = ctL and y = [3L S /L one obtains again the scaling 
form (|l|) with surface stiffness 6=1 and m* = 1. 

Now we turn to the dynamics. For an exact calculation of the transition 
amplitudes (ffi) one may either derive difference equations in t and x as in 



or use the Bethe ansatz [10]. It simplifies matters if (a) we pass to a 
continous-time formulation by taking the limit p — > and t — > oo such that 
t = 2tp remains fixed and (b) by taking periodic boundary conditions (i.e. 
by studying the transfer matrix Tp = T^T instead of T). Neither of these 
simplifications is expected to be of physical relevance to the system || 0|. 
Note, however, that w 2 (x,r) now measures fluctuations in the height dif- 
ferences h(x,r) — h(—L,r). By taking this infinite time limit the evolution 
operator Tp becomes the evolution operator exp (—Ht) of the symmetric 



simple exclusion process ||11|| , but defined on a finite lattice with 2L sites. 
H = — | J2k=-L+i e k is the quantum Hamiltonian of the isotropic ferromag- 
netic Heisenberg model in one dimension. 

In order to derive an exact expression for the height fluctuations ( p!U|) as 
a function of x and r we define the translationally invariant states 

IM = 4= E 44+riO) (l<r<L). (12) 
v £L k= _ L+1 

Note that \r), \ — r) and |r + 2L) are identical and (r\r) = 1 + 5 r> L. Using 
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translational invariance of H one can write with these conventions 



x-l L-l 

w 2 (x,t)=2 !](£ + £ -r)(0,r|exp(-#r)|0,2i?) + 2(a; + L)-(:E + L) 2 

r=-L+l R=l 

(13) 

The matrix elements c R (r, r) = (0, r\ exp (—Ht)\0, R) satisfy the differential- 
difference equation 

-§-c R (r,r) = c R (r+l,r) + c R (r-l,T)-2c R (r,T) (r > 1) 

I ( 14 ) 
£c R (r,r) = c R (r + l,T) - c R (r,T) (r = 1) 

with initial condition c R (r, 0) = <5r,i?(l + &r,l)- We note in passing that one 
may obtain these equations of motion directly from a master equation formu- 
lation of the problem |12] , indicating that the usual stochastic formulation of 
the exclusion process is isomorphic to the continuous-time limit of the model 
defined above. Taking into account the periodicity and reflection properties 
of the states |0,r) one finds 

oo 

Cr{t , r) = e~ 2r Y (-fr-fi+(2L-l)m(2r) + I r+R -l + (2L-l)m 

(2r)) (15) 

m=— oo 

where I p (x) are modified Bessel functions. With (|^) and 

oo i N-1 

£ I p+mN (z) = ^Y eZCOS(27Tk/N)e ~ 2mkp/N ( 16 ) 

we finally obtain 

™ ^ T >~ 2L-1 ^ V ( } costtA;/(2L-1) J 1 - cos27rfc/(2L - 1) +e(X) 

(17) 

where e(x) = (1) for x even (odd). 

It is instructive to compare this exact result with that obtained from a 
phenomenological approach based on a Langevin description [O - the dy- 
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namics of the interface are assumed to be described by the following additive 
noise Langevin equation 



dthj = -(Vi + hj+i - 2h j) + Vj(t) (18) 



where rjj(t) is the usual gaussian white noise. One easily obtains 
°W) = (1 - (-l)>cos2, fa /(2L)) , _ . (19) 



The similarity between the two expressions (|17 ) and ( |1£|) is striking, con- 
sidering the vast simplifications inherent in the phenomenological Langevin 
model. 

It is interesting to take the scaling limit of the above expressions, i.e. we 
take x, t, L — > oo with a = x/L and u = t/L 2 fixed. We then find from both 



(O) and (|19D. 

lim =2 f U ds{ Q^)_ e ^ a _i)/2)) (20) 

where 9^{v) = 9^{v,q) is the theta function with nome q = exp(— tc 2 u). We 
obtain the following asymptotic forms of the scaling function: for arbitrary 
a and u^oowe have 

l im ^M = l(l- tt 2 )- L (l + cos7ra)e- A (21) 



L^oo L 4 



and for a = and aClwe have 



lim — = 2a/m/7t . (22) 



The fact that the scaling function evaluated from the exact dynamics is 
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identical to that obtained from the Langevin description ([Tj^) is very sur- 
prising, since the latter approach completely neglects the strong dynamical 
constraints of the original model, i.e. the restricted possible values of neigh- 
bouring height differences. 

We may gain some insight into this result by considering the correspond- 
ing particle dynamics for the Langevin description (|18|). The appropriate 
particle picture is one of symmetric diffusion with no exclusion, the parti- 
cles to be interpreted as units of height gradient. One may derive an exact 
Langevin equation for this non-exclusive particle process which then may be 
mapped to (]T8|) in the height variables. This is done by Taylor expanding the 
master equation for the distribution P({rij},t) (where rij is the unrestricted 
occupation number at site j), thus deriving a Fokker-Planck equation in the 
Stratonovich representation. The corresponding Langevin equation for the 
occupation numbers may then be mapped to the original height variables 
reproducing (|I8| ). It is of interest to compare this to an exact Langevin equa- 
tion for the exclusion case that we have studied in this paper. Using the 
stochastic Grassmann variables introduced in [12| one may derive an exact 
Langevin equation for the 'height' variables {fj} of the form 

dtfj = \(fj-i + f j+ i ~ 2/0(1 + 2 1 /\(t)) (23) 

(these variables are of Grassmann type, but are simply related to the original 
height variables {hj}). The essential difference between this exact Langevin 



equation and that given in ( 18f) is the appearance of multiplicative noise, 
which appears to be irrelevant in the scaling regime. Therefore the equiva- 
lence between the scaling functions for the restricted and unrestricted inter- 
face models may be thought of as an equivalence between symmetric diffusion 
with or without exclusion - even for non-zero values of the scaling variable 
u. 
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We have yet to derive a full dynamical version of the magnetisation ([[]) 
for our model, as was done in |13| for the Langevin dynamics with additive 
noise. To this end, one has to calculate all higher even moments of the height 
fluctuations. (The odd moments all vanish due to spin reversal symmetry.) 



As regards our model and the phenomenological model of Ref. [ 13] , it is 
unclear whether the identity in the scaling forms of the height fluctuations 
derived above persists beyond the second moments. However, having shown 
that in the scaling limit both the static magnetisation m(x,y) and w 2 (x,t) 
coincide, we conclude that the simplified Langevin dynamics with additive 
noise represent a qualitatively and quantitatively adequate approach to this 
problem in large but finite systems in the scaling region. 

We note that our discussion of the dynamics is limited to the continuous 
time limit defined above. One expects no difference in the scaling form ( |20| ) 
for non-zero hopping rate p (except that the natural scaling variable is now 
t/L 2 as opposed to t/L 2 . However, the continuous time limit defined by 
(1 — p) — > 0,t —>■ oo with t' = t(l —p) fixed has qualitatively distinct scaling 
behaviour since it corresponds to a 'relativistic' limit of the dynamics f9|. 
The behaviour of the interface fluctuations for this case will be discussed 
elsewhere [ H ] . 

In summary, we have investigated a particular stochastic model of a finite 
interface for which the dynamics act on the microscopic degrees of freedom, 
rather than on the phenomenologically specified, coarse-grained microscopic 
variables as is the case in Langevin theories. Nevertheless, results obtained 
so far for our model agree with the earlier Langevin theory, encouraging its 
continued use in appropriate contexts where the truly microscopic theories 
are not analytically amenable. 
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List of Figure Captions 



Fig. 1: The Wulff profile seperating regions of opposite magnetisation. The 
lower (upper) region has negative (positive) magnetisation. 

Fig. 2: The mapping between the restricted interface and the particle ex- 
clusion process. In a) we indicate the initial condition of a flat interface 
corresponding to alternating sites in the particle model being occupied. In 
b) we show a possible interface configuration at some later time and the 
corresponding particle occupancies. The indicated flips in the interface cor- 
respond to particles hopping on the lattice. 

Fig. 3: Allowed vertex configurations in the six- vertex model and their Boltz- 
mann weights. Up-pointing arrows correspond to particles, down-pointing 
arrows represent vacant sites. In the dynamical interpretation of the model 
the Boltzmann weights give the transition probability of the state represented 
by the pair of arrows below the vertex to that above the vertex. 
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